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ABSTRACT 

Multiply-imaged quasars and active galactic nuclei (AGNs) observed in the mid-infrared (MIR) range are commonly assumed to be 
1^ unaffected by the microlensing produced by the stars in their lensing galaxy. In this paper, we investigate the validity domain of this 

assumption. Indeed, that premise disregards microlensing of the accretion disc in the MIR range, and does not account for recent 
progress in our knowledge of the dusty torus which have unveiled relatively compact dust emission. To simulate microlensing, we 
first built a simplified image of the quasar composed of (i) an accretion disc whose size is based on accretion disc theory, and (ii) 
a larger ring-like torus whose radius is guided by interferometric measurements in nearby AGNs. The mock quasars are created in 
the 10'*''^ - 10"** erg/s (unlensed) luminosity range, which is typical of known lensed quasars, and are then microlensed using an 
^ inverse ray-shooting code. We simulated the wavelength and size dependence of microlensing for different Icnscd image types and 

^ ^ fraction of compact objects projected in the lens. This allows us to derive magnification probabilities as a function of wavelength, as 
I well as to calculate the microlensing-induced deformation of the spectral energy distribution (SED) of the lensed images. We find 

O that microlensing variations as large as 0. 1 mag are very common at 1 1 fim (typically rest-frame 4 fim). The main signal comes from 

^ microlensing of the accretion disc, which may be significant even when the fraction of flux from the disc is as small as 5% of the 

C/3 total flux. We also show that the torus of sources with Li,oi <, lO'" erg/s is expected to be noticeably microlensed. Microlensing may 

thus be used to get insight into the rest near-infrared inner structure of AGNs. Finally, we investigate whether microlensing in the 
mid-infrared can alter the so-called i?cusp relation that links the fluxes of the lensed images triplet produced when the source hes close 
^ to a cusp macro-caustic. This relation is commonly used to identify massive (dark-matter) substructures in lensing galaxies. We find 

^ that significant deviations from R^nsp may be expected, which means that microlensing can explain part of the flux ratio problem. 

Key words. Gravitational lensing: micro, strong, quasars: general 
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• 1 . Introduction changes in the flux ratios when a system is observed at a sin- 

. . . . . , 1 1 gle epoch. A flux ratio anomaly is assessed when none of these 

O Gravitational lensmg of quasars is now a well estabhshed tool scenarios account for the observed flux ratios. 

used for a wide variety of astrophysical and cosmological appli- . , i r ..i. • • r ^.i. i ■ ^ 

I ^. , D c J • OAAT 1^ u 1.1 oAA/i ^ popular cxplauatiou for the origm of the anomaly is extra 

^-H cations (Claeskens & Surdei 2002; Kochanek et al. 2004; Treu ^ -c ^- c i j • j * c 

• • ^mrv r. _ 1 '-l^v1^v^ at j i ri i • (de)-magnification ot ouc Icuscd imagc duc to the prescucc ol 

2010; Bartelmann 2010). Mass models of lensing galaxies gen- , ° j, ^ ■ „ ,r.4 ir.8 i . j i 

J ^ J ■ c ^- J ^ -ii- a clump ot mass typically 10 - 10 Mq located close in pro- 

eraUy reproduce the lensed image configuration down to rmlli- . . . . • rf,. . n j n- i • 

■' ^ ,_, 1 <^ 1 N , lection to that image. This ettect, often called milli-lensing 

arcsec accuracy (Chantry et al. 2010; Sluse et al. 2012a), but , . ,,,f , „ „ i • inm n u o 

■' , ^ij a ^ \ .u 1 J • or meso-lensmg (Wambsganss & Paczynski 1992; Baryshev & 

*-j often fail to reproduce the flux ratios between the lensed images „ ir.r.T -k, o o i -j ^nnc,^ ■ a ■ ^ ^ 

d u .1 -irvrvox t j j i v £ lu Ezova 1997; Mao & Schneidcr 1998), IS uowadays au important 

(Yonehara et al. 2008). Indeed, because the light rays from the ^ , „ u- .u . c ■ u . / • i • 

, ,. ,, , . /, , ■ tool for probing the amount of massive substructures in galaxies 

ditterent lensed images cross ditterent regions ol the lensing ^ ^t o -mi a oaai i-. i i e u i oaao c o 

, ^, ..... j-rr . 1 • (Metcalf & Madau 2001; Dalai & Kochanek 2002; Evans & Witt 

Le7;an 1 sTrdei r998"Falco et a" Sndso b^ ^003; Keeton et al. 2003; Vegetti et al. 2010; Metcalf & Amara 

age (Jean & Surdej 1998, Falco et al 1999). Ihey can also be 2012; Xu et al. 2012; Inoue & Takahashi 2012). To use this 
subiect to the microlensing produced by the stars in the lens- . • • 1 j ..1. ..1. -ui 

1 o 1-. r J 1 1 mn T • .11 r,or,^ techniQue, it IS ucccssary to exclude the other possible causes 
ing galaxy (Chang & Refsdal 1979; Irwm et al. 1989). On the r a \ , r,,/ . . . , ^ . t. • r 
.11 J .1 • f- • • i_-i-._ \c.i L- J -.1 of flux perturbations. The mam one is microlensing. It is fun- 
other hand, the mtnnsic vanabmty ol the source, combmed with , c . n- 1 u . 
^, ^. . 1 i_ . .11 . • ... • damentally ol the same nature as nulh-lensing, but it varies on 
the time delay between the lensed images, mtroduces spunous i f .1. . j ^ . 1 . . ■ ^ j 
■' o > t- ^jjjjg scales ol months or years (instead ol at least centuries), and 

* Microlensing maps and source profiles used for our simu- only affects relatively compact sources of projected size similar^ 

lations are available in electronic form at the CDS via anony- 

mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via ' Detectable microlensing can be observed in sources with projected 

http : //cdsweb .u-strasbg . fr/cgi-bin/qcat?J/A+A/ sizes 10 times larger than r]o (Refsdal & Stabell 1991, 1997). 
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to or smaller than one microlens Einstein radius tjo ~ 0.01 pc 
(Sect. 3). For this reason, flux ratio anomalies detected at ra- 
dio wavelengths are considered as the most robust because the 
sources are believed to be larger than 770 (but see Koopmans & 
de Bruyn 2000). Unfortunately, radio flux ratios are measurable 
only for a small sub-sample of lensed quasars bright enough at 
radio wavelengths (about one fifth of the known objects), and 
other proxies of the intrinsic flux ratios are desirable. 

Flux-ratios measured at 8-1 1 yum (corresponding to ~ 2.9- 
4fim rest-frame for typical lensed AGN), are considered as a 
promising alternative to estimate unbiased flux ratios because 
most of the emission is assumed to originate from the large size 
torus (Chiba et al. 2005; Minezaki et al. 2009). Recent progress 
in our knowledge of the dust tori of quasars and active galactic 
nuclei question the generality of this assertion. 

The distance of the dust torus to the central engine of nearby 
AGNs has been obtained using reverberation-mapping, which 
infers the distance to the central black hole from the measure- 
ment of the time-lag between optical and near-infrared (NIR) 
flux variations (Glass 1992; Suganuma et al. 2006, and refer- 
ences therein). Suganuma et al. (2006) compiled the available in- 
ner radii of the dust tori at 2.2 yum for 10 objects, and have found 
that the reverberation mapping radii 7?^^. is consistent with an 
increase with Lp-^, in agreement with the expected dependence 
of the dust sublimation radius with luminosity (Barvainis 1987). 
The reverberation-radius is supposed to be a good approximation 
of the dust sublimation radius, but could be larger depending on 
the exact responsiveness of the dust to the incident continuum 
flux. However, the measured reverberation radius is actually 
smaller by a factor of ~ 3 (Kishimoto et al. 2007) than the naive 
expectation for the sublimation radius using dust properties rep- 
resentative for ISM dust grains (Barvainis 1987). Near-infrared 
(NIR) interferometry has also revealed that the NIR emitting re- 
gion is relatively compact, and that the radius of the torus can 
be as small as the typical Einstein radius of a microlens (Swain 
et al. 2003; Kishimoto et al. 2011a). In addition, we know that 
in the NIR range, there is still a small fraction of the emission 
originating from the accretion disc of the AGN. This means that 
microlensing could be important in the NIR-MIR range. 

Different mass models of the lensing galaxy which reproduce 
the lensed images position, often lead to different predicted flux 
ratios. It is therefore not sufficient to have a flux measurement 
free of extrinsic contamination to identify a flux ratio anomaly. 
Fortunately, for some particular lens-source configurations, the 
flux ratio between lensed images can be predicted in a nearly 
model-independent way. In particular, for the so called "cusp- 
configuration" systems (Fig. 1), the unsigned magnification of 
the saddle-point lensed image^ is expected to be equal to the 
sum of the magnifications of the two (adjacent) minimum im- 
ages. After the seminal work of Mao & Schneider (1998), it 
has been realised that this relation breaks down in the presence 
of massive substructures in the line-of-sight towards one of the 
lensed images, providing a simple technique to identify flux ratio 
anomalies (e.g. Bradac et al. 2002; Keeton et al. 2003). 

In this paper we investigate the efifect of microlensing on the 
NIR-MIR emission of lensed quasars accounting for the refined 
picture we have for the dust torus. We also explore the role of 
the accretion disc emission contributing a measurable fraction of 
the MIR-NIR emission. In Sect. 2, we present our model of the 



~ Lensed images are created at a minima, a maxima or saddle point(s) 
of the arrival time-surface associated the wavefront originating from the 
source. The image parity is calculated based on the magnification tensor 
(e.g. Blandford & Narayan 1986; Schneider et al. 1992). 
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Fig. 1. Example of quadruply-imaged lensed quasar where the lensed 
images are in "cusp configuration" (J 1131-1 231, from Claeskens et al. 
2006). In such a system, the flux of image A (saddle-point of the arrival 
time surface) is expected to be equal to the sum of the fluxes of B and 
C (minima of the arrival time surface). 

source. In Sect. 3, we explain how we proceed to perform the mi- 
crolensing simulations. The results of our simulations are given 
in Sect. 4 and discussed in Sect. 5. In Sect. 6, we investigate 
if the amount of microlensing taking place at MIR wavelengths 
can affect the detection of flux ratio anomalies produced by sub- 
structures, with a particular emphasis on their influence on the so 
called /?cusp relation. In Sect. 7, we study the effect of microlens- 
ing on the SED of the lensed images. Finally we conclude in 
Sect. 8. 

2. Brightness profile of the source 

The continuum emission from AGNs in the range 1-5.5 ;um 
(rest-frame) comes from two components: (i) the emission orig- 
inating from the compact accretion disc, and (ii) the extended 
emission associated with the dust torus. In Sect. 2.1 & 2.2, we 
introduce the relation used to model the brightness profiles of 
these two components. In Sect. 2.3, we explain how we fix the 
relative contribution of each component and how this varies with 
wavelength. 

2. 1 . Accretion disc emission 

We consider that the accretion follows a standard accretion disc 
model with fy oc v'^-' (Shakura & Sunyaev 1973). Kishimoto 
et al. (2008) showed observational evidence of the validity of 
this model in the near-infrared range. We assume the validity 
of this model up to ~ 5 fim (rest-frame). Within this model, the 
half-light radius^ of the disc is given by: 

R^fJA) = 2.37 10'^ cm ( ( ^^t" I (1) 

^ We assume R'^p = 2.44 , with R,i the disc size. This relation is only 
an approximate relation strictly valid if the disc inner radius R^^ « 
R'ljl (See e.g. Abolmasov & Shakura 2012, Sect. 2.3) 
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where /trest is the rest-wavelength of the emission, Mbh is the 
mass of the central black hole, Lboi is the bolometric luminosity 
of the object and tj^s is the accretion efficiency of the disc. In the 
following, we assume a typical accretion rate efficiency /^eff - 
0.1, and a typical Eddington ratio Lboi/^edd = 1/3. Since we fix 
the Eddington ratio, we can derive Mbh fi"om Ledd- Under these 
assumptions, the size of the continuum (accretion disc) emission 

2/3 

is proportional to the bolometric luminosity L^^^^ . 

Some microlensing studies of the X-ray to optical size of the 
accretion disc suggest that: (i) the optical accretion disc is larger 
by one order of magnitude than the accretion disc prediction, and 
(ii) the size of the disc R'^J^ oc A^^^, with /? > 3/4, in disagree- 
ment with the accretion disc theory (Pooley et al. 2007; Morgan 
et al. 2010; Blackburne et al. 2011). However, other microlens- 
ing studies agree well with accretion disc theory (e.g. Bate et al. 
2008; Eigenbrod et al. 2008; Floyd et al. 2009). Possible ex- 
planations to these results involve deviations of some discs from 
the standard disc, under-estimated Eddington ratios, "contami- 
nation" of the optical continuum emission by reprocessed and/or 
scattered light from the disc, systematic errors associated to the 
analyses (Abolmasov & Shakura 2012; Lawrence 2012). How- 
ever, these results concern the optical range, and to date, there 
has been no claim of discrepancy in the (less well explored) NIR- 
MIR range (Agol et al. 2009). 

2.2. Dust torus 

The dust torus of the quasar is modelled as a ring-hke structure. 
For the inner radius of the ring, we adopt the reverberation ra- 
dius in rest-frame K-band (2.2//m) measured by Suganuma et al. 
(2006). Following Equ. 3 of Kishimoto et al. (2007), which 
quantifies the fit of the reverberation radius as a function of L 
done by Suganuma et al. (2006), we can relate the inner radius 
of the torus to the UV luminosity of the AGN: 

/ 6yLv(5500A) \'^' 
^■-^"■"'i lO^^erg/s ) 

We use Lboi = 9.26 v Lv,(5500A) (Shen et al. 2008) to derive the 
bolometric luminosity. 

The other two parameters that define our model are the sur- 
face brightness profile and the outer radius ^out of the torus. The 
surface brightness profile of the torus is poorly known. The- 
oretical considerations support a power-law surface brightness 
(Barvainis 1987) but current data do not allow one to properly 
test this prediction over the NIR (rest) regime. This is however 
not very important for our purposes, because microlensing de- 
pends only weakly on the exact surface brightness as far as the 
same half-light radius of brightness profile is used (Mortonson 
et al. 2005). We therefore use a uniform profile between /?in 
and /?out- Most of the simulations are done for a face-on disc. 
We comment in Sect. 5.3 on the effect of inclination. Interfero- 
metric measurements of local Type 1 AGNs show that the torus 
is relatively compact at 2.2 jim, with ^out/^in ^ 2 (Fig. 7 of 
Kishimoto et al. 2011b). At larger wavelengths (i.e. 8.5 and 13 
fim), the torus tends to be less compact (except for the brightest 
objects) and /?out/^in > 5 is conceivable. Since our observa- 
tional knowledge on the change of the size of the torus both with 
wavelength and luminosity is still quite limited, we choose not 
to parametrise the variation of the size of the torus with wave- 
length. Instead, in this work we consider two extreme scenarios. 
A compact torus with Rom/Rin = 1-2 and a more extended torus 
with /?out/^in = 5. Under these assiunptions, half of the light of 



the ring is between /?i„ and RfJ!^^ ~ 1.1 /?in for the compact torus 
model and ~ 3.6 /?in for the extended model. 

2.3. Compound model 

In order to build a realistic model of the quasar source, it is 
necessary to estimate what fraction of the accretion disc con- 
tributes to the total flux at a given wavelength. The modelhng 
of the SED of the nearby AGNs used for interferometric stud- 
ies suggests that the continuum contributes 14-29% of the to- 
tal emission at 2.2 jum (Kishimoto et al. 2011a). Therefore, we 
assume a fiducial contribution of the accretion disc to the flux 
of /ad = 0.2 at 2.2 jum. In order to derive /ad at other wave- 
lengths, we use /v oc v'^^ for the accretion disc and we consider 
a sum of two black-body emissions for the dust. The use of 
two black-body components (a "high" and a "low" temperature 
component) for the dust emission provides a reasonably good 
description of the near-infrared SED of Type 1 AGNs. We use 
a "high" temperature component BBl with a fiducial tempera- 
ture T - HOOK, and a "low" temperature component BB2 with 
T = 300 K, which dominates the dust emission at /I > 7 fim. The 
relative contribution of each component is given by the ratio be- 
tween the efi'ective area of the corresponding emitting region. In- 
terferomettic measurements by Kishimoto et al. (201 la) indicate 
that the effective area of the BB2 region emission is 100-1000 
times larger than the BBl region emission (Table 7 of Kishi- 
moto et al. 2011a). We choose a typical ratio of efl'ective area 
^2/^^! = 400 leading to the following parametrisation of the 
torus flux F'"""' = 400Bv(v,r = 300 K) + By{v,T = 1400 K). 
Our model of the spectral energy distribution of the source is 
shown in Fig. 2. 

In the following, we investigate the effect of microlensing at 
rest-frame wavelengths of 1, 2.2, and 4.4 jim. We choose these 
wavelengths because they roughly correspond to the observed 
A'-band, L-band and 1 1 jum of the lower redshift lensed quasars 
(1.1 ^ z ^ 1.5). At these wavelengths, the contribution of the 
cold T = 300 K component of the dust remains less then a few 
percent. It is therefore unnecessary to model explicitly its sur- 
face brightness separately. We only include it in the calculation 
of fAoi-i), to account for the fact that the accretion disc becomes 
more and more diluted with increasing wavelength. 

Most of the known lensed quasars have unlensed bolometric 
luminosities Lhoi in the range [10"*^ - 10^^'^] erg/s. Therefore, 
we decided to study microlensing for 3 different luminosities: 
Liow = 1044-2 erg/s, L,„ed = 1045 erg/s and Lhigh = lO^*- erg/s. We 
show in Table 1 the characteristic size of the source, in units of 
the Einstein radius t/q (Eq. 3), for these luminosities. 

3. Microlensing 

The scale-length at which the source is typically affected by mi- 
crolensing is the Einstein radius of the microlenses, given by: 

where the D are angular diameter distances, and the indices o, 
1, s refer to observer, lens and source, respectively, and (M) 
is the aver age mass of the microlenses. The value t]o ~ 3 
xlO^^cm V{^)/0.3 Mq ~ 0.01 pc is not derived for precise 
source and lens redshifts but is roughly the average 770 obtained 
for known lensed quasars (Mosquera & Kochanek 201 1). All 
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Fig. 2. Spectral Energy Distribution (SED) of our source model be- 
tween 0.5 and 11 yum (rest frame). The accretion disc (AD) emission 
is shown with a solid red line, the "hot" component of the torus with a 
dashed green line, the "cold" component with dashed-dotted blue line 
and the total SED with thick black line. 

Table 1. Sizes of the accretion disc (AD) and of the dust torus used in 
our model for three different luminosities (col. #1). Col. #2 displays the 
half-light radius of the accretion disc (in units of the Einstein radius rjo 
for (M> = 0.3 Mq). The next columns characterise the properties of the 
torus. Col. #3, 4, 5 give the type of torus, its inner radius R,„, and the 
width of the ring. These quantities are given for A^^^t = 2.2/jm. Only 
R'^J^ oc A*'^ varies with the wavelength (Eq. 1). 



tended sources as a function of position. The spatial resolution 
we considered fixes the lower bound of our luminosity range, 
Lbo\ = lO^"*'^ erg/s. For that source, we have R1/2 = 0.02 770 at 
2.2 yum. Our lower resolution microlensing maps are needed to 
be able to estimate the amount of microlensing of the extended 
torus. These maps are however still too small to calculate mi- 
crolensing for the extended torus when L = 10"'^erg/s. In that 
case, Rom - Rm > 100 t/o- This size of the torus is too large to be 
affected by significant microlensing and we therefore assume no 
microlensing of that model of the torus. 

The amplitude and the probability of (de)-magnification of 
a source because of microlensing depends on the image type, 
or more precisely, on the convergence k and shear y at the lo- 
cation of the lensed images. Microlensing also depends on the 
fraction of matter in smooth and compact form at the projected 
location of the lensed images. Since we cannot sensibly ex- 
plore the whole range of plausible values of k and y, we have 
decided to focus on values of k and y typical of saddle-point 
and minimum image observed in cusp-like systems. Our fidu- 
cial "saddle-point" image has (a:, y) - (0.47, 0.57), and our fidu- 
cial minimum has (/f, 7) = (0.421,0.5), coiTesponding to an 
average macro-magnification |yL(.„,d| - 19.8 and |yL(,„,„| - 11.8. 
These values match those obtained for image A and B of the 
lens Jl 131-1231 (Fig. 1) using an elliptical Singular Isothermal 
Ellipsoid + external shear (SIE+7) to model the lensing galaxy 
(Sluse et al. 2012a). Recent microlensing studies of X-ray flux 
ratios in lensed quasars indicate that the most likely fraction of 
compact matter in lensing galaxies is = k^/k - 0.07, at the 
location of the lensed images (Pooley et al. 2012). This is the 
fiducial value we use in this paper. Because there is evidence 
for a larger fraction of compact objects in some systems, we also 
test for Kt - 0.3 and Kt = 1. 



Lboi (erg/s) Rfjl (t]q) torus type Rj^ (tjq) R^^ - Rj,, (t/q) 
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along this work, we use this typical value to convert linear source 
sizes (derived in cm) into Einstein radii. 

Microlensing of the source is estimated in the following 
way'*. First we build a micro-magnification map using the 
inverse-ray-shooting code developed by Wambsganss (1990, 
1999, 2001). This map gives, for a pixel-sized source, the micro- 
magnification as a function of the source position. Since we 
want to study microlensing affecting a source with both a com- 
pact (continuum) and an extended (torus) structure, we need 
maps with large extent and relatively high resolution. In prac- 
tice, we use maps of IOO770XIOO/70, with a spatial resolution of 
0.01 770/pixel for Lboi < 10"*5 erg/s and maps of 25077ox250/7o, 
with a spatial resolution of 0.05 Tyo/pixel otherwise. Second, fol- 
lowing the methodology described in Sect. 2, we build the sur- 
face brightness of the source on a grid as fine as the microlens- 
ing map. Third, we convolve our image of the source with 
the microlensing map. Each pixel of that convolved map pro- 
vides the amount of micro-magnification associated to the ex- 

The microlensing maps and source profiles used in this work are 
made available via CDS. 



4. Results 

The amount of microlensing affecting a macro-lensed image of 
flux F2 is quantified by the factor fj., which accounts for the extra- 
magnification induced by microlensing such that Ff''^ - juF^. 
Hereafter, following e.g. Wambsganss (1992), we calculate for 
various source and lens properties, the probability density func- 
tion (PDF) for the quantity Am = 2.5 \ogQi). With that con- 
vention. Am > (/i > 1) corresponds to a magnification of the 
source and Am < is a demagnification. 

For a fixed lensed image property (corresponding to a given 
combination of k, y), three sets of parameters influence the prob- 
ability of microlensing in the NIR-MIR range: (a) the physi- 
cal characteristics of the source, namely its luminosity, and the 
compactness of the dust torus; (b) the relative fraction of com- 
pact and smooth dark matter in front of the lensed image; (c) the 
variable contribution with wavelength of flux originating from 
the torus and from the disc. 

4.1. Changing the source properties 

We first investigate how the luminosity of the source and the 
compactness of the torus modify the probability of microlens- 
ing. We follow the prescriptions of Sect. 2 for /I = 2.2 yum 
to calculate the brightness profile of the source. In order to 
gain intuition on the relative effect of the dust and of the ac- 
cretion disc, we show in Fig. 3 the PDF of microlensing for 
each component of the emission separately. The three left-most 
panels show the distribution for the accretion disc, for a com- 
pact torus and for an extended toiTis, while the two right-most 
panels show the distribution for the compound source models. 
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The accretion disc contributes to 20% of the compound source. 
In each panel, the PDF for three different source luminosities 
log(Lboi/erg/s) = [46,45,44.2] are compared. Panel (a) shows 
the PDF for a saddle-point image, and panel (b) for a minimum 
image. In both cases, we fix the fraction of compact matter to 
Kt = 0.07. We compare panels (a) and (b) in Sect. 4.2 while here 
we focus on the results of the saddle-point image (i.e. panel (a)). 

First, this figure confirms that the compact emission from 
the accretion disc can be significantly microlensed for all lumi- 
nosities considered. Second, we see that if the torus is as com- 
pact as suggested by near-infrared interferometric data, the prob- 
ability for it to be significantly microlensed is non-negligible, 
even for intrinsically bright sources. On the other hand, if the 
torus emission is more extended, the probability of getting the 
torus microlensed by > 0.1 mag drops to almost zero. Finally, 
we see from the rightmost panels that the small contribution of 
the flux from the accretion disc governs the probability distribu- 
tion of microlensing, such that independently of the torus model, 
there is always a significant probability for the source to be mi- 
crolensed in that wavelength range. 

Another interesting feature unveiled by that figure is the exis- 
tence of a "cut-off" in the PDF of the compound model for large 

demagnification (the PDF drops sharply at Am 0.25 mag). 

This effect is discussed in Sect. 5.1. 

4.2. Changing the image type 

In this section we compare the PDF of microlensing for a saddle- 
point image (Fig. 3a) and a minimum (Fig. 3b). We can see 
that the accretion disc (left-most panel) emission in a saddle 
point image has a much larger probabiUty to be strongly de- 
magnified (Am < -1.1 mag) than the minimum image. This 
difference in behaviour between saddle-point and minimum is 
a well known effect and we encourage the interested reader to 
consult Schechter & Wambsganss (2002) and Saha & Williams 
(2011) for a deeper understanding of this effect. Interestingly, 
we see that this difference becomes weaker with the increasing 
source size. This is particularly clear for the torus, when com- 
paring the PDF for the saddle-point and for the minimum image. 
This similarity of the PDF between the minimum and the saddle- 
point persists for the compound model. The origin of this effect 
is the "suppression" of the large demagnification in saddle-point 
images. Indeed, when strong demagnification occurs, the flux of 
the accretion disc gets diluted in the flux of the torus to such an 
extent that only a weak demagnification can be observed. It is 
however important to note that the suppression of large demag- 
nification in the near-infrared is only valid in the context of mi- 
crolensing by a population of solar-mass-like microlenses. If the 
Einstein radius of a microlens is larger than the torus size, i.e. the 
situation encountered in case of milli-lensing by substructures in 
the lensing galaxy with M > lO^Mg, then large demagnifica- 
tion of the torus (and therefore of the compound model) can still 
occur. 

4.3. Changing the lens properties 

The fraction of compact matter in front of a strongly lensed im- 
age is in the range of 5 to 10% (Mediavilla et al. 2009; Pooley 
et al. 2012), but seems to be several tens of percent in other sys- 
tems or even reach 100% for the Einstein Cross (Dai et al. 2010; 
Poindexter & Kochanek 2010). In Fig. 4, we show the PDF for 
microlensing of a compound source at A = 2.2 fim for three dif- 
ferent values of /c*. The results are displayed for three different 



source luminosities, and for the two torus models. It can be seen 
that a decrease of k» leads to a shift of the peak of the PDF to- 
wards significant demagnification. This is again a well known 
effect for saddle-point images (Schechter & Wambsganss 2002). 
However, it is important to notice that this effect is mostly caused 
by the sizable contribution of the accretion disc to the total flux 
(see Sect. 5). 

4.4. Changing the wavelength and /ad 

The appearance of the source varies significantly from 1 to 
4.4 yum. The main change comes from the decreasing contribu- 
tion of the accretion disc emission to the total flux with increas- 
ing wavelength. According to our model, the decrease remains 
however relatively small, and 10% of the flux still comes from 
the disc at 4.4jum (see Sect. 5.1). On the other hand, the ac- 
cretion disc size also increases with wavelength. We display in 
Fig. 5 the PDF for three wavelengths and luminosities. We see 
that the shape of the PDF changes with increasing wavelength 
(peaking at larger A) because of the progressive dilution of the 
flux of the accretion disc in the flux of the torus, not because of 
the increase in the size of the torus and of the disc with A. 

Since we see that the accretion disc is the main cause of mi- 
crolensing in the rest-NIR, we investigate how the probabihty 
of microlensing changes with /ad- For this purpose we fix the 
accretion disc and torus sizes to their values at 2.2 yum, and calcu- 
late the probability of observing microlensing (de)magnification 
\Am\ > 0.1 mag for different values of /ad- Figure 6 shows the 
results for different values of k^. We see that the probabihty to 
observe \Am\ > 0.1 mag is already significant for /ad ~ 0.05. 
The compactness of the torus only weakly changes the results, 
except when Lboi < 10^^ erg/s and when /ad is small. This is in- 
deed the regime where the compact torus can be significantly mi- 
crolensed such that the probabihty of microlensing gets boosted 
despite the small contribution of flux from the accretion disc. 

4.5. Typical probabilities of microlensing 

Table 2 sunnmarises our results. We provide the probability for 

a source to be microlensed by more than 0.1 mag for differ- 
ent source luminosities and fraction of compact objects /<:, in the 
lens. We give those probabilities for the accretion disc (AD) at 
2.2 /zm, for our two different kinds of tori (i.e. compact and ex- 
tended torus, see Table 1), and for the compound source model 
(AD+torus) at 3 different rest-frame wavelengths. 

Table 2 shows that microlensing of the dust torus amounts 
more than 0. 1 mag for the compact torus model and sources with 
Lboi ^ 10"^^ erg/s. Microlensing of the torus remains below 0.1 
mag for the extended torus. The probability of significant mi- 
crolensing of the (compact) torus increases with a:,, conversely 
to the behaviour observed for the accretion disc. These results 
are quahtatively independent of the image type (i.e. minimum or 
saddle-point). 

The probability of microlensing for the compound model is 
governed by the relative contribution of the accretion disc flux 
to the total flux. As shown in Fig. 6, values of /ad ~ 0.05 are 
sufficient to enable microlensing Am > 0.1 mag. As soon as 
the accretion disc contribution is large enough, the probability 
of microlensing depends on both the image type and /c*. We see 
at e.g. 4.4 jum that saddle-point images show a larger probabihty 
of microlensing Am > 0.1 mag for small Kt, while minimum 
images show a large probability for Kt - I. 
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Fig. 3. Probability density function (PDF) for microlensing Am at 2.2 yum for the accretion disc (AD), the torus (compact and extended model) 
and compound model (torus+accretion disc). The distributions are given for three different luminosities log(L(,oi/erg/s)=44.2, 45, 46. Note that in 
our convention Am > corresponds to magnification. 




Fig. 4. PDF for microlensing Am of a saddle point image at 2.2 fim for three different fractions of compact matter (k,) and for three different 
luminosities log(Lboi/erg/s)=46, 45, 44.2. Solid lines correspond to the compact torus and dashed lines to an extended torus. Note that in our 
convention Am > corresponds to magnification. 
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Fig. 5. PDF for microlensing Am of a saddle point image at three different wavelengths A = l.,2.2,4.4yum, and three different luminosities 
log(L(,oi/erg/s)=46, 45, 44.2. Solid lines are for a compact torus and dashed lines for an extended torus. Note that in our convention Am > 
corresponds to magnification. 
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Fig. 6. Probability to observe microlensing (de)magnification |Am| > 0. 1 mag as a function of the fraction of the flux /ad originating from the 
accretion disc. The solid line is for a compact model for the torus and the dashed line for an extended model. The blue, red and green curves 
correspond to different source bolometric luminosities. Panels from left to right correspond to k, = 0.07, k, = 0.3, k, = 1.0. Probabilities are 
calculated for A^^^t = 2.2 yum. 

Table 2. Probability of observing microlensing (de)magnification | A»i| > 0. 1 mag in a saddle-point (minimum) image, for three source luminosities 
L, and three values of k,. Probabilities are calculated for a pure accretion disc (AD) at 2.2 /im, a pure torus (compact or extended), and for 
compound source models as they would be observed at A = 2.2 yum, 1.0 //m, 4.4 yum. 



Saddle-point image 



L = 10'*'' erg/s 


L = 10*^ erg/s 


L = lO"*"*-^ erg/s 


0.07 0.3 1.0 


0.07 0.*3 1.0 


0.07 0.3 1.0 


AD 0.91 0.89 0.84 


0.95 0.91 0.85 


0.95 0.91 0.86 


Compact 0.00 0.00 0.00 
1.0 yum 0.94 0.89 0.83 
2.2 yum 0.59 0.49 0.30 
4.4 yum 0.06 0.12 0.07 


0.03 0.21 0.23 
0.94 0.90 0.84 
0.74 0.60 0.41 
0.48 0.41 0.27 


0.27 0.50 0.48 
0.94 0.91 0.84 
0.74 0.64 0.55 
0.54 0.53 0.49 


Extended 0.00 0.00 0.00 
1.0 yum 0.94 0.89 0.83 
2.2 yum 0.59 0.47 0.31 
4.4 yum 0.05 0.08 0.04 


0.00 0.00 0.00 
0.94 0.90 0.84 
0.74 0.57 0.36 
0.46 0.27 0.08 


0.00 0.00 0.00 
0.94 0.91 0.84 
0.75 0.61 0.39 
0.56 0.36 0.17 


Minimum image 


AD 0.84 0.87 0.86 


0.89 0.89 0.88 


0.89 0.90 0.88 


Compact 0.00 0.00 0.02 
1.0 yum 0.86 0.88 0.86 
2.2 yum 0.23 0.42 0.39 
4.4 yum 0.00 0.08 0.14 


0.01 0.25 0.35 
0.87 0.89 0.87 
0.35 0.53 0.54 
0.13 0.37 0.43 


0.17 0.48 0.57 
0.87 0.89 0.87 
0.41 0.60 0.59 
0.27 0.51 0.57 


Extended 0.00 0.00 0.00 
1.0 yum 0.86 0.88 0.86 
2.2 yum 0.21 0.40 0.37 
4.4 yum 0.00 0.04 0.06 


0.00 0.00 0.00 
0.87 0.89 0.87 
0.30 0.51 0.43 
0.06 0.17 0.12 


0.00 0.01 0.04 
0.87 0.89 0.87 
0.33 0.53 0.47 
0.08 0.30 0.28 



5. Discussion 

The non-negligible probability of observing a microlensing sig- 
nal of at least 0.1 mag in the rest-frame range 1-4.4 yum, im- 
plies that microlensing cannot be neglected in the study of lensed 
quasars flux ratios from A'-band up to 1 1 yum (observed frame). 
We discuss hereafter how to link our simulations to observed flux 
ratios in lensed quasars. 

5. 1 . Contribution of tlie accretion disc 

The fraction of the flux originating from the accretion disc (/ad) 
controls the amount of microlensing in the NIR-MIR. In our 
model, we used a fiducial value /ad = 0.2 at /Irest = 2.2 yum. 
This value is typical of what is observed in local AGNs (i.e. /ad 
in the range 14-29%). We show in Fig. 7 how /ad varies with 
wavelength. It is important to realise that fAoi^) may change 



depending of the real slope of the power law accretion disc emis- 
sion in the near-infrared. In our model, we used a standard index 
ffv - 1/3. This choice is supported by spectropolarimetry of 
a sample of local AGNs (Kishimoto et al. 2008). It contrasts 
with the redder spectral slope ay 0.4 observed in the optical- 
ultraviolet domain (e.g. Vanden Berk et al. 2001). Such a redder 
slope, if valid at larger wavelengths, would lead to an even larger 
contribution from the disc in the NIR-MIR and would increase 
the contribution of microlensing in that range. We illustrate the 
effect on fAoi'^) for different values of a,, in Fig. 7. 

The various PDF shown in Figs. 3, 4, 5 display a sharp cutoff 
for demagnification but not for magnification. The origin of this 
cut-off is caused by the effective increase of contrast between 
the flux from the accretion disc and from the torus in case of 
large demagnification. When this situation occurs, the effective 
fraction of the flux coming from the accretion disc drops sig- 
nificantly, such that the total flux gets dominated by emission 
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Fig. 7. Variation of the fraction of the flux coming from the accretion 
disc (/ad) with rest-frame wavelength. /adW is shown for three dif- 
ferent power law index Oy of the accretion disc emission. For the three 
curves we have fixed /AD(2.2//m) = 0.2, the fiducial value used in this 
work. 

from the torus and that extreme demagnification has no effect 
anymore. In case of micro-magnification, the inverse occurs, 
and the total flux becomes larger than if the torus was the only 
source of emission. 

The cut-off value of the PDF can be derived by expressing 
the microlensing of the compound model: 

Am = 2.5 log [/ad /^ad + /^toms (1 - /ad)] • (4) 

From this equation, we see that when the accretion disc is signifi- 
cantly de-magnified, ^uad — > 0, we have Am 2.5 log[/itorus (1 - 
/ad)]- Therefore, when the demagnification of the torus remains 
small (/itoms — > 1), the minimum demagnification of the source 
wiU be Am = 2.5 log(l - /ad)- On the other hand, there is no 
formal upper limit to the maximum magnification, and the latter 
is only governed by the maximum amount of microlensing of the 
accretion disc (see Refsdal & Stabell 1991, 1997), and by /ad- 

5.2. Extended vs compact torus 

We find that noticeable microlensing of the torus is likely to take 
place with amplitudes as large as 0.15-0.2 mag for objects with 
Lboi ~ lO^ -^ erg/s if a compact model for the torus is assumed. 
Based on interferometric studies of nearby AGNs, this model is 
appropriate for rest-frame wavelengths ~ 2.2 jum. At a rest-frame 
wavelength of 5.5 yum or larger^, the hot component of the torus 
gets several times more extended than at 2.2 fim, especially for 
Lboi ^ lO'*^ erg/s (see Fig. 7 of Kishimoto et al. 2011b). There- 
fore, the simulation considering an extended torus may be more 
appropriate when deaUng with A ~ 4.4 /im (rest-frame). For 
larger luminosities (Z-boi ~ lO''^ erg/s), the torus remains com- 
pact (/?out/^in ~ 2), but the half-light radius is significantly larger 
than a microlens Einstein radius such that microlensing does not 
exceed 0.05 mag (Fig. 3). Consequently, it is correct to consider 
that at 4.4 /im (rest-frame), most of the microlensing originates 
from the disc. This is however with the assumption that there is 
no evolution with time of the torus properties. This is supported 

' There is unfortimately no measurement in the rest-frame range 2.2- 
8.5yum. 



by a recent study of the properties of the 1-10/zm SED of AGNs 
detected in the COSMOS field which finds no strong evolution 
of the SED with redshift (Hao et al. 2012). 

5.3. The effect of inclination 

Type 1 AGNs are not always face-on as considered in our model 
and inclinations of the polar-axis with respect to our line-of-sight 
as large as 60 deg are plausible. Inclination will reduce the ob- 
served area covered by the source and thus increase the prob- 
ability of important microlensing. We show in Fig. 8 the PDF 
for a saddle-point image at /Irest = 2.2 jim, and an inclination of 
the polar axis of 60 deg. We see that the microlensing distribu- 
tion of the torus gets broadened by a factor of about two. On 
the other hand the probability of large (de)-magnification of the 
accretion disc gets slightly increased, especially for the intrinsi- 
cally brighter sources. These two effects lead to a broadening 
of the PDF of the compound model and to a shift of the peak of 
PDF towards demagnification for Lboi = 10^^ erg/s. 

We have also examined the effect of inclination for a min- 
imum image and k» = 0.07. In that case, the distribution of 
magnification of the accretion disc remains roughly unchanged, 
but the magnification distribution of the torus gets broader, simi- 
larly to the situation encountered for the saddle-point. However, 
this effect remains too small to significantly modify the PDF of 
the compound model associated to the minimum. 

5.4. Existing observations and simulations 

To our knowledge, there are only a few observational studies of 
the SED of lensed quasars in the NIR-MIR range. Agol et al. 
(2000) were the first to study the MIR ffux ratios in the Ein- 
stein cross = Q2237-(-0305. They aimed to disentangle the syn- 
chrotron and dust emission origins of infrared AGN flux. Their 
analysis demonstrated that the MIR emission (observed frame) 
in Q2237-I-0305 originates from a region larger than one mi- 
crolensing Einstein radius r]o (for (M) = 0.3 Mg), ruUng out 
pure synchrotron emission. This result was later confirmed by 
Wyithe et al. (2002) who explored a larger family of macro- 
lens models. Recently, Agol et al. (2009) used updated MIR 
flux ratios obtained with the IRAC camera on Spitzer to improve 
the study of the SED of Q2237-H0305. Their observations (see 
their Fig. 7) showed a clear chromaticity of the flux ratios be- 
tween 0.6 and 11 /im (observed frame). A similar trend was ob- 
served by Minezaki et al. (2009) using ground-based MIR-data. 
Agol et al. (2009) demonstrated that in order to reproduce the 
chromaticity of the flux ratios, the emission in the rest-range I- 
4.2 //m had to originate from two regions, a hot torus and an 
accretion disc with a temperature profile compatible with pre- 
diction from the standard thin accretion disc model. Two other 
quadruply lensed quasars have been observed in the MIR by 
Chiba et al. (2005). These authors, searching for milli-lensing 
in these systems, also suggest significant compact emission in 
the system B 1422+231, but not in PGl 115+080. The doubly- 
lensed quasar HEl 104-1805 has also been observed in the MIR 
and its flux ratios in the context of microlensing have been anal- 
ysed by Poindexter et al. (2008). The lack of chromaticity of 
the flux ratios in that system suggests that the emission down 
to ~ 2.2 jim was dominated by an extended emission possibly 
originating from the torus. The data however do not rule out 
that a minor fraction of that emission comes from the disc. Ob- 
servations of lensed systems above 2.2 /im are generally sparse. 
Chromaticity of the flux ratios in the NIR (up to 3.8 jum observed 
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Fig. 8. Effect of inclination on the probability density function (PDF) of microlensing of a saddle-point image. The dashed lines show the PDF of 
Fig. 3 (face-one quasar) and the solid lines show the corresponding PDF when the polar axis of the quasar is inclined by 60 deg. The distributions 
are given for log(Lboi/erg/s)=44.2, 45, 46. Note that in our convention Am > corresponds to magnification. 



frame) has been investigated in six lensed quasars by Fadely & 
Keeton (2011). Some systems show a "plateau" in the flux ra- 
tios between 2.2 and 3.8 /im while others show clear diff'erences 
likely caused by microlensing or milli-lensing (Fadely & Keeton 
2012; Sluse et al. 2012b). 

A recent paper of Stalevski et al. (2012b), which appeared 
after the start of this work, investigated the generic possibility 
of microlensing of the dust torus. Our work differs in several 
respects from Stalevski et al. (2012b) and complements their re- 
sults. The main conceptual difference between the two studies 
relates to the model of the source. The work from Stalevski et al. 
(2012b) focuses on the torus while we also address the effect of 
the accretion disc on the flux ratio. On the other hand, they use 
a torus model derived from 3-D radiative transfer simulations of 
a two-phase medium (Stalevski et al. 2012a; Baes et al. 2011, 
2003). Our model is simpler, but it is based on observational 
measurements of the torus size. Owing to the use of Rin equal to 
Barvainis sublimation radius (Barvainis 1987), they effectively 
use a torus with an inner radius 3 times larger than observed 
in nearby objects (Suganuma et al. 2006), and than our model, 
for the same L. For these reasons, direct comparison of results 
have to be careful. Because of the difference of the definition of 
inner radius, their source with L - 10'^ Lq = 3.9 x 10'*''erg/s 
should match our source with Lboi ~ lO^^erg/s. Contrary to 
Stalevski et al. (2012b), we do not have a model accounting for 
the increase of Rout with wavelength but we consider two ex- 
treme models for the torus, a compact and an extended torus. At 
A = 2.2 jum, their model compares to our compact torus model, 
and at A - 4.4 //m, to our extended torus model. From their 
Fig. 4 (left panel), for a lens^ at z = 0.5, a maximum microlens- 
ing magnification of 0.06 mag at 2.2 fim is expected, in agree- 
ment with our Fig. 3. We cannot perform the comparison at 
4.4 yum because we assumed no microlensing of the extended 
torus for the high luminosity source. 

5.5. Variations witin time 

Contrary to milli-lensing, the amplitude of microlensing is ex- 
pected to vary on a human time scale (Schechter & Wambs- 
ganss 2002). This time scale depends on the relative trans- 
verse velocity between the source, the microlens and the ob- 
server (Kochanek 2004). Typical source plane transverse veloc- 
ities are v, ~ O.OST/o/yri and can reach v, ~ 0.15?7o/yr in the 

^ This choice of redshift is needed to roughly match the Einstein radius 
with our fiducial value of rjo. 



most favourable cases (Mosquera & Kochanek 201 1). We show 
in Fig. 9 an example of an expected microlensing light curve in 
a saddle-point image for a source with Lboi - lO'*'*'^ erg/s, and a 
compact torus. The figure shows the light curve separately for an 
accretion disc alone (AD), a torus model alone, and for a more 
realistic composite model (see Sect. 5.1). The light curves are 
displayed for three different wavelengths: l.l/jm, 2.2/im, and 
4.4 yum. It can be seen that the most drastic variations with wave- 
length are generally associated to caustic crossing-like events. In 
such cases, chromatic variations will be observed from the ultra- 
violet, where the accretion disc dominates, up to the near/mid- 
infrared range. If the source is located farther away from a caus- 
tic, then the chromaticity will be observed nearly exclusively in 
the near/mid-infrared range. In that case, the chromatic change 
is mostly caused by the variable fraction of flux from the ac- 
cretion disc to the total flux. Another interesting feature is the 
suppression of large demagnifications when the flux of the torus 
takes over the flux of the accretion disc (cf Eq. 4). Because the 
microlensing of the accretion disc does not show important chro- 
maticity in this situation, we could in principle use large demag- 
nification events to measure /ad(^)- 

As outlined previously, the microlensing signal from the 
torus is relatively weak. In the most favourable case of our com- 
pact torus model, flux variations might be observed with am- 
plitude as large as 0.1 mag on time scales of typically years to 
decade(s). For a brighter source and/or for more extended torus, 
the microlensing of the torus rarely reaches more than 0.05 mag 
and time scales for significant variations reach several decades. 

6. Consequences for flux ratio anomalies 

We investigate in this section the strength of the anomalies pro- 
duced by microlensing in the MIR. To quantify this effect, we 
simulate microlensed flux ratios for minimum and saddle point 
images, and calculate the quantity /?cusp- This quantity is com- 
monly used to identify flux ratio anomalies caused by substruc- 
tures (e.g Mao & Schneider 1998; Keeton et al. 2003; Xu et al. 
2012). The magnified fluxes Ia,b.c of the bright saddle-point im- 
age A and two minima images B and C in a cusp-configuration 
(Fig. 1) can be used to calculate: 

„ _ |1-/b//a-/c//aI 

1+/b//a + /c//a- ^ ^ 

In the ideal case, Rcusp ~ 0. Significant deviations from zero 
or from the value estimated based on a simple macro model is 
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Fig. 9. Example of microlensing light curve for a saddle-point image 
(with K, = 0.07). The microlensing-induced flux variations (in mag) as 
a function of time are displayed for a source with L(,oi = lO*'*-^ erg/s, 
and a compact torus model. The microlensing event duration scales 
with the inverse of the net transverse velocity v, expressed in Einstein 
radi/year. The upper panel shows a portion of the microlensing map 
from where the light curve is extracted. Blue (dark) region correspond 
to demagnification and green (light) region to magnification. The other 
panels show the light curve for the accretion disc (AD), for the torus, 
and for the composite source (AD+torus). The results are presented for 
three different wavelengths, 1 fjm (light blue), 2.2 fim (light green), and 
4.4 fim (red). Note that the torus model is the same for each wavelength. 

considered as a good sign of an anomaly. For our fiducial values 
of (a:, y) of the saddle-point and minimum images, we calculate 
^cusp - 0.068. We show in Fig. 10 the distribution of values of 
^cusp in presence of microlensing, for different source luminosi- 
ties, and for fAoi^ - 4.4 fim) - 0.2, 0.1 and 0.02. We see that a 
large deviation from the expected value of Rcusp may be observed 
in the MIR range, especially for an intrinsically fainter system. 

7. Consequences on the SED 

We show in Fig. 1 1 (left panel) how the SED of a source with 
Lbo\ ~ 10"*^ erg/s can be modified if the accretion disc is (de)- 
magnified due to microlensing by typically Am ~ 0.6 mag at 
A ~ 1 yum. We have picked two particular events, but the situation 
remains very similar if random events are selected. Although 
this amount of microlensing is relatively large, it is not uncom- 
mon and is indeed observed in several lensed systems (Sluse 
et al. 2012b). In this example, we have assumed an extended 
torus and our fiducial value of /ad- Because of the assumed 
source luminosity and extended torus model, only the accretion 
disc is significantly microlensed. We see that in case of demag- 
nification, the peak of the the SED at ~2.2 fim gets more pro- 
nounced because of the increased contrast between the torus and 
the disc, while in case of magnification, the contrast is decreased 
and the bump in the SED gets nearly suppressed. For an even 
larger magnification, this bump can disappear completely. The 
right panel of Fig. 1 1 shows the ratio between the microlensed 
and unlensed SED. This is qualitatively similar to the chromatic 
changes which would be observed in the flux ratio between 2 
lensed images, one being much more aff'ected by microlensing 
thant the other. It is ubiquitous that the main change of the SED 
takes place below /lrest=2 /vm. The changes of the SED depend 



on fAoi/i) but also to some extent on the characteristics of the 
microlensing event. This precludes the direct use of ratios of 
SEDs of different lensed images to infer /ad- The time varia- 
tion of the event might help to characterise the microlensing but 
the monitoring time scales needed are discouragingly large (see 
Sect. 5.5). 

The main interest of studying a microlensed rather than an 
unlensed quasar relies in the change of contrast between the 
torus and the disc emission induced by microlensing. The most 
useful targets in this context are quadruply imaged systems. In 
such a case, microlensing is independent in each lensed image, 
and it is likely, that at a given time, one of the lensed image will 
be affected by large microlensing while another one will only 
be weakly microlensed. The estimate of the amplitude of mi- 
crolensing of the accretion disc in the images can be based on 
spectroscopy, comparing the flux ratios in the continuum and in 
the broad emission lines of the different pairs of lensed images 
(Sluse et al. 2012b; Sluse et al. 2011; Hutsemekers et al. 2010). 
Alternatively, a comparison of the SED of the different lensed 
images could also be informative since similar SEDs suggest 
low microlensing (de)magnification. The most microlensed im- 
age, in case of significant demagnification of the accretion disc 
emission, could be used to study the emission from the torus. 
In that case, microlensing acts as a natural coronograph, hid- 
ing a fraction of the disc emission. On the other hand, in the 
case of significant magnification of the accretion disc emission, 
study of the intrinsic shape of the disc will be possible up to 
longer wavelengths. Finally in a situation of monitoring of a 
micro-magnification event, comparison of the observed signal 
with microlensing simulations could be used to derive the in- 
trinsic disc temperature profile up to rest-NIR wavelengths (e.g. 
Anguita et al. 2008; Eigenbrod et al. 2008; Blackburne et al. 
2011). 



8. Conclusions 

We have studied microlensing of large sources taking place in 
images of strongly-lensed quasars when observed in the near- 
/mid-infrared wavelength range (typically from 2 to 1 1 /vm). Our 
simulation considers objects with intrinsic bolometric luminosi- 
ties^ in the range lO"*'*'^ - lO"^*" erg/s. Our two-component model 
of the quasar emission allows us to distinguish between mi- 
crolensing of the dust torus, whose size scales as L**^, and of 
the accretion disc, whose size scales as L^^^ (for a fixed Edding- 
ton ratio). We find that microlensing with amplitude > 0. 1 mag 
is expected to be common in the considered wavelength range. 
Most of the expected signal is due to microlensing of the accre- 
tion disc. The latter cannot be neglected as soon as the disc con- 
tributes to more than 5% of the total flux, a situation expected 
to be common. On the other hand, microlensing of the toiTis 
becomes noticeable for the lowest luminosity quasars we con- 
sidered, and at rest-wavelengths typically < 2.2 fj.m, i.e. where 
the torus is still compact and emission is dominated by the hot 
(T ~ 1400 K) dust component. We may expect targets with 
Lbo\ ^ 10"^"* erg/s to show even stronger microlensing of the torus, 
offering the perspective of studying the latter with microlensing. 
Such faint systems are however uncommon among the known 
lensed quasars. This might change with upcoming all sky sur- 

' It is important to bear in mind that the accretion disc source size de- 
pends also on the black hole mass. Because we have assumed a constant 
Eddington ratio L/Ledd = 1/3, we could express our results as a function 
of the source luminosity only. Therefore, our results indirectly depend 
on the Eddington ratio. 
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Fig. 10. Probability density function i?cusp at 4.4 yum (rest-frame) for /ad = 0.2 (left), /ad = 0.1 (middle), /ad = 0.02 (right). The distributions 
are given for three different luminosities log(L(,oi/erg/s)=44.2, 45, 46. The upper panel shows the distribution for a compact torus and the bottom 
panel for an extended torus. The value of R^usp in absence of microlensing is R^usp = 0.068 and is depicted with a vertical dashed line. 



veys like LSST or through dedicated lensed-quasars hunting pro- 
grammes^. 

It is well known that the probability distribution of mi- 
crolensing of a compact source is different for images which are 
saddle-points or minima of the arrival-time surface (Schechter 
& Wambsganss 2002; Saha & Williams 201 1), the former being 
more prone to demagnification than the latter From the analysis 
of microlensing of the torus, we could see that the differences 
are less significant for more extended sources. More interest- 
ingly, the magnification distribution of the compound source is 
very different from the one obtained for the disc alone. The main 
characteristic of the PDF of the compound model is the suppres- 
sion of large demagnifications. The reason of this suppression 
is simply that, when the disc is demagnified, its flux gets diluted 
in the total flux to such extent that only a weak demagnification 
is observed. If /ad is the fraction of flux from the disc, and if 
the torus is assumed to be unlensed, the maximum demagnifica- 
tion of the compound model is given by Am = 2.5 log(l - /ad)- 
Conversely, in case of magnification, the contrast between the 
torus and the disc decreases, and large amplitude microlensing 
remains possible. 

One of our aims was to estimate if microlensing could ex- 
plain flux ratio anomalies in the NIR-MIR. For that purpose, 
we have calculated the expected distribution of i?cusp for quasars 
of different luminosities, assuming a cusp system similar to 
J1131-1231 (Sluse et al. 2003). We found that for almost all 

Any such programme, if based on colour selection criteria, will have 
to account for colour changes induced by microlensing to be efficient. 



the situations considered, Rcusp may get significantly perturbed 
due to microlensing. A possibility to identify if an anomaly is 
caused by substructure or not might be through the study of the 
spectral energy distribution (SED) of the lensed images up to 
4yum (rest-frame). We have shown that in case of important mi- 
crolensing, the SED could be significantly modified. This sug- 
gests that the comparison of the SED of the different lensed im- 
ages should allow one to disentangle between microlensing and 
milli-lensing. In case of milli-lensing by substructures of mass 
of at least IQ^Mq, the dust torus should be significantly magni- 
fied (because the Einstein radius of the substructure is similar 
to the size of the near-infrared torus), leading to less severe de- 
formations of the SED. An estimate of the differential effects 
between the accretion disc and the torus (and within the torus) 
in the case of milli-lensing is needed to assess the absence of 
change of the SED in that case. 
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